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ABSTRACT 

We study the large-scale anisotropic two-point correlation function using 46,760 luminous red galax- 
ies at redshifts 0.16 - 0.47 from the Sloan Digital Sky Survey. We measure the correlation function 
as a function of separations parallel and perpendicular to the line of sight in order to take account of 
anisotropy of the large-scale structure in redshift space. We find a slight signal of baryonic features in 
the anisotropic correlation function, i.e., a "baryon ridge" corresponding to a baryon acoustic peak in 
the spherically averaged correlation function which has already been reported using the same sample. 
The baryon ridge has primarily a spherical structure with a known radius in comoving coordinates. 
It enables us to divide the redshift distortion effects into dynamical and geometrical components and 
provides further constraints on cosmological parameters, including the dark energy equation-of-state 
With an assumption of a flat A cosmology, we find the best-fit values of fim 



0.218to;o37 and 



Qb = 0.047lg;Q;^g (68% CL) when we use the overall shape of the anisotropic correlation function of 
40 < s < 200 h~^ Mpc including a scale of baryon acoustic oscillations. When an additional assump- 



tion of Qbh — 0.024 is adopted, we obtain fJoE = 0.770 



+0051 n _ 

0.040 ana w — 



-0.93 
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These constraints 



are estimated only from our data of the anisotropic correlation function, and they agree quite well 
with values both from the cosmic microwave background (CMB) anisotropics and from other comple- 
mentary statistics using the LRG sample. With the CMB prior from the 3 year WMAP results, we 
give stronger constraints on those parameters. 

Subject headings: cosmological parameters — cosmology: observations — galaxies: distances and 
redshifts — large-scale structure of universe — methods: statistical 
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1. INTRODUCTION 

Recently, baryon acoustic oscillations have been ob- 
served in the large-scale structure of the universe. These 
observations include an analysis of the two-point corre- 
lation function (2PCF) of the Sloan Digital Sky Survey 
(SD SS) Luminous Red Gala xy (LRG) spectroscopic sam- 
ple (jEisenstein et al.l 120051 E05, hereafter), the power 
spectrum of the Two-Degree Field (2dF) Galaxy Red- 
shift S ur vev (iCole et al.ll2005h a nd the SDSS LRG (lliltsil 
I2006airbl : iTegmark et al.l l2006l : iPercival et all l2007allbD . 
and the angular power spectrum of the SDSS LRG 
samp l e with photometri c redshifts (jPadmanabhan et al.l 
120071: iBlakeet al.l 120061 ). There is a hi nt of acoustic 
oscill ations in the SDSS quasar sample (jYahata et al.l 
l2005f ). These analyses have established the ability of 
the baryon acoustic oscillations to constrain cosmolog- 
ical parameters competitively and c omplementarily with 
the C MB fe.g..lSpergelet alll2007f) and Type la super- 
novae (jRiess et al.lll998l: iPerlmutter et al.lll999D . 

These previous analyses of the baryon oscillations, 
however, use angle-averaged 2PCFs, £,{s), or angle- 
averaged power spectra, P{k), where s and k are the 
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separation and wavenumber in redshift space, respec- 
tively. A certain amount of information is lost when 
anisotropi cs of structure are i g nored . In their pioneer- 
ing work, lAlcock fc Paczvhskil ()1979D proposed that ge- 
ometrical anisotropics in redshift space of the high-z 
universe can be used as a probe of t he cosmological 
consta nt. iMatsubara &: Sutol (|1996D and iBallinger et aT] 
(|1996l ) pointed out that the anisotropy of the 2PCF 
and the power spectrum in redshift surveys can con- 
strain the dark energy components. Recently, methods 
which directly use anisotropy of the baryon acoustic os- 
cillations have be en theo retically developed for both the 
power spectrum (iHu fc Haiman 2003; Sco & Eiscnstcip 
120031: iGlazebrook fc Blake 2005; Sco & Eisenstcin 200"^ 
and the 2PCF (Matsubara 2004 1. These approaches use 
functions of two variables: separations parallel and per- 
pendicular to the line of sight. Observationally, esti- 
mations of such two-variable functions are noisier than 
one-variable functions. The baryon acoustic signature 
in the large-scale structure is weak even in one- variable 
statistics when presently available samples of galaxies 
are used (see, e.g., E05). Therefore, methods that di- 
rectly treat the anisotropy of the bar yon acoustic fea- 
ture require very lar ge surv ey v ol ume (Eisenstcin et alj 
119991: IMatsubara fc Szalav 200ll) . [Peacock et al. C200lf 
and iHawkins et al.l ( 20031 ) m easured the 2PCF with two 
variables from the 2dFGRS (jColless et al.ll2001f) and de- 
tected the detailed signature of large-scale coherent in- 
fall, and as a result were able to constrain a value of 
(3 ~ fl'^/b which parameterizes linear redshift distor- 
tions. However, they used information from scales much 
smaller than the acoustic scale because of limited survey 
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volume and did not analyze the effect of geometrical dis- 
tortion. On the other hand, the correlation analyses of 
the 2dF QSO survey placed constraints on the cosmo- 
logical constant (Hoyle et al. 2002; Outram et al. 2004; 
da Angela et al. 2005; see also Ross et al. 2007). Their 
analyses still focused on smaller scales than the baryon 
acoustic scale. 

In this paper we analyze the anisotropic 2PCF, includ- 
ing the baryon acoustic peak in the large-scale structure 
and constrain relevant cosmological parameters. We use 
a spectroscopic sample of the SDSS LRG, which is the 
most useful sample for our purpose. Our analysis differs 
from previous studies of baryon acoustic oscillations in 
the LRG s ample in that, ow ing to a theoretical devel- 
opment bv iMatsubaral (|2004f ). we take into consideration 
the fully two-dimensional feature in the 2PCF to detect a 
geometrical distortion effect; this is the first cosmological 
application of the two-dimensional acoustic peaks. 

Before proceeding to the next section, we take 
note of nonlinearity on large scales and the scale de- 
pendence of the galaxy biasing. The nonlinear ef- 
fects on the baryonic features that appear around 
100 h~^ Mpc play an important role in taking account 
of percent-level cosmo l ogy (e.g., [M eiksin et al. 119991 : 
ISeo fc EisensteinI 120051 : iJeong fc KMnatsul l2006f) . Re- 
cent work also suggests that the scale-dependent biasing 
poses a serious probl e m in analyzing galaxy surveys (e.g., 
Blanton et aT '2006': 'P ercival et aT '2007 bl: ISmith et all 
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However, for the sake of simplicity, in this paper we con- 
sider only the large-scale clustering and assume the bi- 
asing to be scale-independent and linear. 

The structure of this paper is as follows. In § [2] we 
describe the SDSS LRG sample used in our analysis. We 
then, in § [21 measure the anisotropic 2PCF in redshift 
space and estimate its covariance matrix. In § |4] we out- 
line the cosmological parameter dependence on the mod- 
eled 2PCF, including dynamical and geometrical distor- 
tions. Cosmological parameters are constrained by the 
measured anisotropic 2PCF in §[5j In §[B]our conclusions 
are given. 

2. THE SDSS LRG SAMPLE 



The SDSS (lYork et al.l [2000t iStouehton et~al] [200l 
is an ongoing imaging and redshift survey which uses 
a dedicated 2.5m telescope, a mosaic CCD camera, 
and two fiber-fed double spectrographs (|Fukugita et al.l 
Il996t iGunn et al.l I l998. 20 0i). After image process- 
ing ("Lupton et al."2001; Stoughton et al." 2002t|Pier et all 
[2003; Ivezic ct al. 2004; Tucker et al. 200l) and calibra- 
tion (Hog g et all I2001t ISmith et all I2003 ). the spectro- 
scopic targets of LRGs are selected from the imaging data 
accord ing to the algorithm described bv lEisenstein et aTl 
(I2001D. The tiling al gorithm for the fibers is found in 
lBlantonet"aII(|2003a[ ). 

For our analysis, we use 46,760 LRGs over 3853 deg 
in the redshift range from 0.16 to 0.47, which is the same 
sample as the one used in previous analyses of SDSS 
LRG clustering (Zehavi et al. 2005; EO S). The sky cov- 
erage is the same as lss_samplel4 (Blan ton et al.ll2005f ) 
and is similar t o that of the p ublicly available SDSS 
Data Release 3 (jAbazaiian et al.i.2004 ). The galaxies in 
the sample have rest-frame ^-band absolute magnitudes 




Fig. 1. — Illustration of geometrical quantities used in this paper. 

-23.2 < Mg < -21.2 {Ho = 100 km s"^ Mpc^^) with 
K + E corrections of passively evo lved galaxies to a fidu- 
cial redshift of 0.3 (jBlanton et al.i.2003b. ). The comoving 
number density of the sample is close to constant out 
to 2; = 0.36 (i.e., volume Hmited) because of the narrow 
absolute magnitude cut, and drops thereaft er due to the 
flux limits (see Fig. 1 of lZehavi et~alll2005( l. The radial 
and angular selection functions, fiber collisions, and un- 
ob served plates are mo deled using the method described 
in IZehavi et al.l (|2005D : E05 also provides the details of 
the sample. 

3. ANISOTROPIC 2PCF OF LRGS 

3.1. Measuring the LRG 2PCF 

The 2PCF is measured by comparing the actual galaxy 
distribution to a catalogue of randomly distributed 
points in the same r egion, accordi ng to the selection func- 
tion of the survey ([Peebles! 1 1980f ). We count the galaxy 
pairs in bins of comoving separation along and across 
the line of sight, s\\ and s±_ respectively, to estimate the 
anisotropic 2PCF. Our notations of geometric quantities 
are illustrated in Figure [H First, the comoving distances 
to every galaxy, x{z)^ are calculated by assuming a fiat 
universe with fim = 0.3 and I^a = 0.7, where rim is 
the mass density parameter and f2A is the cosmological 
constant parameter. This fiat universe is used only as 
a mapping between the observed space and our analysis 
space: our theoretical modeling also takes this mapping 
into account. The purpose of the mapping is simply to 
avoid having to perform the analysis in strongly distorted 
redshift space. For each galaxy with redshift zi we define 
the separations s\\ and s±^ of other galaxies with redshift 
Z2 according to Figure [1] 

(1) 

(2) 

where is the apparent angle between the two galaxies 
from the observer. This definition of the line of sight is 
not as standard as that of (vp, tt), such as was defined 
by Davis & Peebles (1983; see also Fisher et al. 1994). 

Which galaxy of the pair is chosen as galaxy 1 is arbi- 
trary in our definition. Both galaxies are considered as 



S|| —x{z2) COS 9 — a;(zi), 
s± —x{z2) shi9, 
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galaxy 1, and the line of sight is simply defined as the 
direction toward the galaxy 1. Therefore we count a pair 
of galaxies twice. As a result, there appears a strong 
correlation between certain bins in {s±, sii) space; in 
particular, the bins of opposite sign of S|| are strongly 
correlated. Those bins contain almost identical pairs in 
a small-angle case, 6* <C 1. However, the angle 6 in our 
sample is not always this small, so that those bins con- 
tain different sets of pairs. Still the correlations between 
those bins are strong, which is properly taken into ac- 
count in our parameter estimation below. 

We compute th e anisotropic 2PCF u sing the Landy- 
Szalay estimator (jLandv fc: Szalavl[l993f ). 
, DD-2DR + RR 

^(■'^^'\0 = RR ' (3) 

where DD, RR, and DR are the normalized counts 
of galaxy-galaxy, random-random, and galaxy-random 
pairs, respectively, in a particular bin in the space of 
(s_L, S||). The random catalogue contains about 30 times 
as many points as the real data, and the random points 
are distributed according to the radial and angular se- 
lection functions. The space of (sj_, sh) is divided into 
rectangular cells with As±, As\\ = 10 h^^ Mpc. Each 
galaxy of redshift z is weighted by 1/[1 + n{z)Pw], where 
n{z) is the comoving number density and Pm is the power 
spectrum at a typical scale (jFeldman et al.|[r994[ ). We 
adopt Pw = 40, 000 h~^ Mpc , which is evaluated at the 
baryon wiggle scale and is the same value as E05. We 
have also t ried another va l ue, P yj = 30, 000/i~^Mpc^, 
adopted bv iTegmark et all (|2006D . We found that the 
va lue does not have a str ong effect on the result, as noted 
in lPercival et al.l ()2007bD . because the comoving number 
density of our sample is close to constant at almost all 
scales. 

The resulting redshift-space 2PCF for the observed 
LRGs is shown in the right half of the plane in Figure [H 
The value of ^(s_l, S||) are given by contour lines. There 
is an indication of the baryon ridges of radius about 
100 h~^ Mpc, which is a counterpart of the baryon peak 
in the one-dimensional 2PCF, although the signal is not 
so strong. The anisotropy of the clustering is obvious in 
this figure. When the separation is along the line of sight 
[s± « 0), the clustering is elongated due to nonhnear ve- 
locity dispersions of galaxies. On the other hand, the 
large-scale clustering is squashed along the line of sight 
due t o coherent infalls toward over-dense regions (JKaiserl 
Il987( ). The latter effect is often called the Kaiser's ef- 
fect. A corresponding theoretical pre diction based on a 
linear perturbation theory derived by iMatsubaral (|2000l . 
12004) is shown in the left half of the plane in Figure O 
Although the measured 2PCF is noisy, the linear theory 
can account for the behavior of the 2PCF on large scales. 
The nonlinear velocity distortions are not described by 
linear dynamics, which should be removed in our linear 
analysis below. The detailed comparison clearly needs 
statistical treatment, which we explain below. When the 
anisotropic 2PCF is averaged over the angle, the one- 
dimensional 2PCF ^(s) is obtained. 

3.2. Covariance Matrix 

Because there are strong correlations between different 
bins of the anisotropic 2PCF, it is necessary for statis- 
tically proper analyses to estimate a covariance matrix. 
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Fig. 2. — Contour plots of the redshift-space 2PCF measured 
from the SDSS LRG sam ple (right) and the corresponding ana- 
lytical formula derived by IMatsubaral 1120041 ^ using a linear per- 
turbation theory (left). The dashed black lines show ^ < —0.01 
increasing logarithmically with 0.25 and —0.01 < C < linearly 
with 0.0025. The solid thin lines colored red show < ^ < 0.01 
increasing linearly with 0.0025 and the solid thick ones colored red 
£, > 0.01 logarithmically with 0.25. The baryonic feature slightly 
appears as ridge structures around the scale s = (s^ + s?.)^'^ ~ 

100 h~^ Mpc, and the dashed circle traces the peaks of the baryon 
ridges . For the theoretical predictions, we adopt the best-fit values 
assuming a fiat cosmology, Qm = 0.218, Qt = 0.0473, h = 0.702, 
as = 0.660, b = 1.55, while the fiducial values, rig = 1 and w = —1 
are fixed. We also set the redshift at the origin to be 0.34, which 
is typical in our LRG sample. 

For this purpose, jackknife resampling or bootstrap re- 
sampling (e.g., see Lupton 1993) is often adopted. In a 
cosmological context, however, cosmic variance plays a 
critical role in estimates of cosmological parameters. It 
is uncertain whether these methods can provide a reli- 
able estimator of the cosmic variance because they rely 
only on one observed sample. We first tried to use the 
jackknife method and found that this approach actually 
underestimates the covariance at all the scales, which 
leads to the underestimation of the error bars for cosmo- 
logical parameters when we compare them with a more 
reliable met hod explained be l ow. A similar tendency is 
also seen bv iPope fc Szapudil (|2007t ). 

One of the best ways to estimate the covariance matrix 
including cosmic variance is to use N-body simulations 
to generate many mock catalogs from which the covari- 
ance matrix is calculated. It is necessary to generate a 
larger number of mock samples than the number of data 
points of the statistics to be employed; otherwise, one 
would improperly obtain a singular covariance matrix. 
However, because our analysis of the anisotropic 2PCF 
deals with several hundred data points, it is computa- 
tionally too expensive to produce a sufficient number of 
independent realizations in our case. 

In our analysis we adopt an alternative method, us- 
ing a public seco nd-order La grangian perturbation the- 
ory (2LPT) code (jCrocce et al. 2006). As input ingredi- 
ents, we adopt flm = 0.27, Qa = 0.73, h = 0.7, cts = 0.8, 
256^ particles in a cubic box of side 1600 h^^ Mpc, 
and a transfer function calculated by CMBfast code 
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Fig. 3. — Comparison of the 2PCFs times s^ between the ob- 
served LRGs and the mock catalogs. The horizontal axis is loga- 
rithmic while the vertical axis is linear. The red points show the 
angle-averaged 2PCF of the LRGs and the error bars are from the 
mock catalogs. The dashed line shows the 2PCF averaged over 
the mock catalogs. The baryon peak detected in these plots is ob- 
tained from the integration of baryon ridges in Fig. |2]over angular 
orientation. ,. . 

(|Seliak fc Zaldarriagal[l99l with fit, = 0.045, where h is 
the Hubble parameter normahzed by 100 km s~ Mpc~^ , 
ug, is the rms. of the fluctuations smoothed with a top- 
hat window function of radius R — 8 h~^ Mpc, and fib is 
the baryon density p arameter. To impla nt the galaxy bi- 
asing of the LRGs (jZehavi et al.ll2005l ). we empirically 
select particles with probability proportional to e""'™, 
where Sm is the mass density fluctuation at a position 
of each particle calculated by the University of Washing- 
ton HPCC's public SMOOTH code.^ We choose a = 1.5 
so as to match the correlation amplitude of the observed 
LRGs at scales larger than 40 h~^ Mpc. Then we trim 
the mock catalogs from simulation boxes so as to have 
the same number density and the same survey volume 
as the actual LRG sample. Finally, we generate 2500 
mock catalogs with independent initial conditions, com- 
pute the 2PCF for each, and obtain a covariance matrix 

by 



1 



N 



(C),, ^ Cov(e„ e,) - ^^T^ E(^' - ^"»)(^i - ^^■)' (4) 



1=1 



where N = 2, 500, f ' represents the value of the 2PCF 
of ith bin in Ith realization, and ^t is the mean value of 
£,1 over realizations. The average of 2PCFs from each 
mock catalog agrees well with the observation within the 
1 a errors for both one-dimensional (Fig. [3]) and two- 
dimensional analyses (Fig. [4]) . These mock catalogs are 
constructed solely to estimate errors of the measured 
2PCF of LRGs, so the averaged 2PCF over the cata- 
logs is not used for the following analyses. Increasing a 
changes the amplitude at the scales around the baryon 
peak to match the observation, while it causes more dis- 
crepancy at the small scales. This tendency is, however, 
consistent with the theoretical prediction (see Figure 3 
in E05). 

^ At|http://www-hpcc. astro. washington.edu/tools/smooth. html I 
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Fig. 4. — Anisotropic 2PCFs as functions of two variables, sep- 
arations perpendicular and parallel to the line of sight. The right 
side shows the LRG 2PCF, which is the same as the right one of 
Fig. [2] The left side shows the corresponding averaged 2PCF of 
our mock catalogs. The difference between 2PCFs in each mock 
and their average is used for construction of the covariance matrix. 
The 2PCF from our mock catalogs does not have large deviation 
from that of the observed 2PCFs, even for the quadrupole compo- 
nents. 

We test the full covariance matrix obtained by the 
2LPT method. First, we randomly choose the 2PCFs of 
30 reahzations out of 2500. We regard each of 30 2PCFs 
as an observed LRG sample, calculate the x^ statistics 
by the method described in ij5l constrain the input cos- 
mological parameters, flm and h, and check how many 
realizations contain the input values at 68% and 95% con- 
fidence level. If the amplitude of covariance is reasonable, 
68% of all the contours of 68% confidence levels should 
contain the input parameters. The reason for using only 
30 realizations is that comparing against all the mock 
catalogs is computationally very expensive; it requires 
that x^ be calculated in seven dimensional parameter 
space 2500 times. According to our statistics, anisotropic 
2PCF ^(sj^,S||), there are 22 and 28 realizations which 
contain the inputs at 68% and 95% confidence levels for 
rim in the total 30 realizations, while there are 17 and 28 
for h. We thus conclude that our method of estimating 
the covariance matrix is reasonable, and can be reliably 
applied for parameter estimation in fJSl Figure [5] shows 
the result of the test; we choose to present only 15 real- 
izations because displaying all 30 results makes the figure 
unclear. 

4. THEORETICAL PREDICTIONS 

As a theoretical prediction for the anisotropic 2PCF, 
we a dopt an analytical formula of Matsubara (200(J, 
l2004f ) derived in a ge neral situation taking into account 
the wide-angle effec t (ISzalav et all [l998i) and the high- 
z distortion effect (jMatsubara &: Sutol Il996l ) in linear 
perturbation theo ry. The necessary formula is given in 
iMatsubaral (|2004[ ). 

The mean redshift of the LRG sample is about 0.34 
and the clustering scale which we probe ranges up to 
200 h~^ Mpc. The maximum angle between two points 
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Realization 

Fig. 5. — Display of the reliability of the recovery of 51^ and 
h from 15 realizations. The horizontal axis shows the realization 
number, while the vertical axis shows the value of Om and h and 
the horizontal dashed lines show their input parameters. Solid 
and dashed error bars shows the 68% and 95% confidence levels, 
respectively. Among these 15 realizations, there are 10 and 11 
realizations including input values of f2m and h, respectively, at 
68% confidence intervals , and 14 and 15 at 95% intervals. 

from the observer {0 in Fig. [T|) is ~ 12°. The distant- 
observer approximation is not so accurate at > 10°, and 
therefore the general formula give above, which accu- 
rately includes the wide-angle effect, is preferable for a 
precise analysis of the LRG sample. 

The left side of Figure [2] shows a prediction in linear 
theory of the two-dimensional 2PCF in rcdshift space 
with the central redshift of zi = 0.34. We adopt a 
flat cosmology with f7„ = 0.218, f)fc = 0.0473, /i = 
0.702,(78 — 0.660, h = 1.55, Ug = 1, and w = —1, where 
w = Pde/pde is the equation-of-state parameter for the 
dark energy component and b is the linear bias param- 
eter. The first five values are our best-fit values for the 
two-dimensional 2PCF of LRGs, as described in the fol- 
lowing section, while the last two are the fiducial values. 

Throughout this paper we assume a flat cosmology for 
simplicity. There are seven cosmological parameters in 
our modeling: Vim, h, fif,, Us, w, cts, and b. For the de- 
tails of the dependence of the anisotropic 2P CF on cos- 
mological parameters, see iMatsubaral ()2004[ ). In short, 
there are three kinds of physical effects. The first one is 
the shape of the underlying mass power spectrum, which 
is determined by the f2m, ^6, ^, and Us. The second 
one is the dynamical distortion effect which is generated 
by peculiar velocities of galaxies. Linear, coherent ve- 
locities squash the apparent clustering along the line of 
sight, while nonlinear, random velocities smear the clus- 
tering along the same direction (JKaiseilllQSTt iHamiltonl 
Il992f ). The linear squashing effect depends on the so- 
called redshift distortion factor, (i{z) — f{z)/b{z), where 
f{z) = d\nD/d\na is the logarithmic derivative of the 
linear growth rate D{z) at redshift z, a ^ {l + z)~^ is the 
scale factor, and b{z) is the linear bias factor at redshift 
z. The growth factor depends on flm and w. Since we 
assume a flat cosmology, the density parameter of dark 
energy is given by Ode = 1 — ^m- However, the pa- 
rameter dependence on the growth factor is not so useful 
in parameter estimation, because the overall amplitude 
of the power spectrum characterized by erg is a free pa- 
rameter. Since the parameter dependence of nonlinear 



velocity effect is not analytically given, we do not use 
the nonlinear regime in the 2PCF. The third effect is 
the geometric distortion, which depends on the Hubble 
parameter H{z), and angular diameter distance, Da{z), 
and thus depends on Qm, and w. The dependence of the 
geometric distortion on h vanishes in redshift surveys in 
which distances are measured in units of h~^Mpc. The 
geometric distortion is useful for constraining the dark- 
energy parameters, J^de and w. The baryon ridges are 
isotropic in comoving space and their anisotropy is pri- 
marily due to geometric distortion. 

Finally, we comment on the galaxy biasing and the 
evolutionary effect. Because we consider only the linear 
regime, we assume the biasing to be scale-independent 
and linear. Although the bias parameter b is completely 
degenerate with as in the ordinary one-dimensional 
2PCF, the two-dimensional 2PCF is able to solve this 
degeneracy through the measurement of the redshift dis- 
tortion parameter, f3, which depends on b. We treat b and 
ag, as independent parameters in the following analysis. 
One could choose /3 or berg as free parameters (which are 
more closely related to the measurements), instead of b 
or cTg. The choice of the independent parameters does 
not affect the following result. 

In this paper we consider the measured 2PCF as a 
representative of the function at a mean redshift zi = 
0.34. We therefore simply neglect the effects of evolution 
on clustering and biasing within the sample. Strictly 
speaking, evolutionary effects are not negligible in very 
large redshift surveys which have broad range of redshift 
(e.g.,iYamamoto & Suto 1999). For example, the evolu- 
tion has a significan t effect in the SDSS quasar sample 
(jYahata et al.ll200"5h . In our LRG sample, however, the 
redshift range is relatively small and the signal-to-noise 
ratio of the measured 2PCF is not very high. Indeed, 
the evolutionary effect on the growth factor is about less 
than 20% from the survey edge to the mean redshift, but 
the effects on the anisotropic 2PCF and cosmological pa- 
rameters are negligibly small compared to error levels. 

5. CONSTRAINTS ON COSMOLOGICAL MODELS 

5.1. Setup 

In this section we describe methods and results of 
constraining cosmological parameters by the anisotropic 
2PCF of the LRG sample. We measure the goodness of 
fit, which shows how well assumed cosmological parame- 
ters fit a set of observational data, and the measurements 
are given by the x^ statistics, taking into account the full 
covariance matrix. As described in the previous section, 
we adopt seven free parameters rim (= 1 ~ ^de), ^6, h, 
Hs, w, (Tg, and b, assuming flatness of the universe. 

In comparing observational data with theory, we first 
compute the theoretical 2PCF as a function of {zi,Z2, 9) 
with a given set of parameters. Next we use the fiducial 
parameters of Vl^ = 0.3 and Oa = 0.7 to convert red- 
shifts into comoving distances, which are the same val- 
ues assumed for measuring the distances of each galaxy in 
§ 13.11 Therefore, we compare the theoretical 2PCF with 
the observation in the same comoving space and these 
fiducial parameters do not bias our results of parameter 
estimations. 

The theoretical formula based on the linear perturba- 
tion theory does not reproduce nonlinear gravitational 
effects and nonlinear velocity distortions such as finger- 
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of-God effects. We therefore discard the observed 2PCF 
at scales less than 40 h~^ Mpc. We also do not use the 
data along the hne-of-sight, namely, s± < 10 h~^ Mpc, 
becaus e the line-of-sigh t components of the 2PCF are 
noisy teerns teinl Il994l ) and furthermore deviat e from 
Kaiser's formula even on large scales (:Scoccimarral200^ . 
Finally we perform the analysis for the scale range of 
40 < s < 200 h~^ Mpc. We also adopt a more conserva- 
tive range 60 < s < 160 h~^ Mpc to check the systematic 
effects beyond the linear theory. The numbers of bins in 
2PCF are 574 for 40 < s < 200 h-^ Mpc and 330 for 
60 < s < 160 h^^ Mpc. x^ statistics are then calculated 



as 



X 



'■ie) = J2^^^io)ic-%A^,i9), 



(5) 



where 6 is a set of cosmological parameters to be con- 
strained, A^i{0) denotes the difference between the ob- 
served and theoretical 2PCFs in ith bin, and the sum 
is over the number of bins. The most likely values 
for the cosmological parameters minimize the equation 
([S]). Finally, the likelihood function for the cosmologi- 
cal parameters, C, is proportional to exp(— x^/2) with 
an appropriate normalization factor. Then, for example, 
the 68% confidence interval becomes the region where 
J C dd — 0.68 in the parameter space. 

5.2. ag-b degeneracy 

First, we consider the behavior of the parameters re- 
lated to the clustering amplitude, ag and b. Figure [S] 
plots their joint likelihood functions with contours rep- 
resenting 68%, 95% and 99% confidence levels, where 
w — —1 is fixed and we marginalize over the other four 
parameters flm, ^b, 



h, and ns. We find as = 0.66;^;2i6 



and b = 1.55_q'75 (68% confidence level) for the fit to 
40 < s < 200 h~^ Mpc after also marginalizing them 
over each other. As described in §|4l b and cg are strongly 
coupled as the amplitude of the 2PCF is proportional 
to a product ba-g. The degeneracy is somewhat alle- 
viated from anisotropy of the 2PCF due to dynamical 
distortions which are dependent on /3. It is still diffi- 
cult, however, to independently constrain these two pa- 
rameters without relying on other observations such as 
CMB, higher order correlation analysis, etc. In this pa- 
per we mainly focus on parameter constraints only from 
the 2PCF and consider the joint analysis with the Wilkin- 
son Microwave Anisotropy Probe ( WMAP ) results only 
in §[531 

Therefore, we always marginalize over both as and b 
in the following likelihood analysis. This marginalization 
corresponds to mainly using the shape information in the 
2PCF and discarding the amplitude information. 

We also note that the physical origins of the parame- 
ters b and as are not fully understood. The value of the 
bias parameter b depends on unknown details concerning 
the formation of LRGs, and the value of as depends on 
unknown details about the generation of density fluctu- 
ations in the primordial universe. There is not any re- 
liable theory which robustly predicts the values of these 
parameters. 

5.3. Main results 

We next focus on the four fundamental cosmological 
parameters J7m, flf,, h, and Ug after marginalizing over 
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Fig. 6. — Likelihood contours for the parameters related to the 
clustering amplitude, trg and b. The diagonal panels represent the 
likelihood functions for the individual parameters, where Cm, O;,, 
h, and ris are marginalized over. The bottom left panel represents 
a two-parameter constraint, and each ellipse shows the 68%, 95%, 
and 99% confidence levels from inward. The solid and dashed (col- 
ored red and black) contours are the fits for 40 < s < 200 h~^ Mpc 
and 60 < s < 160 h~^ Mpc, respectively. The best-fit parameters 
are trg = 0.66 and b = 1.55 for 40 < s < 200 h'^ Mpc. 

as and b. Figure [7] illustrates contour plots of the joint 
likelihood functions of two parameters among the four, 
where w is fixed at —1. 
The fits to 40 < s < 200 h'^ Mpc give n 



0.218 



+0.047 
-0.037' 



^b = 0.04731^; 



+0.0157 !, _ n 7r)o+0.187 - 

0160' " — ^-'^^-0.117' ^^^ 

Us = 1.122to;l83 (68% CL). The best-fit values and the 
errors of all the cosmological parameters are listed in Ta- 
ble [TJ The best-fit values of all of the four parameters 
for 40 < s < 200 and 60 < s < 160 h'^ Mpc are con- 
sistent within the range of the 68% error. This result 
suggests that the systematic effects beyond the linear 
theory are small. The accuracy of the constraints on 
Qm and fif, increases by nearly a factor of 2 while the 
uncertainty of h is improved only marginally when the 
wider dynamical range is used. We obtain a relatively 
strong constraint on Qm without fixing the value of fif, 
because the constraint from the linear squashing effect in 
the anisotropic 2PCF is stronger than in the spherically 
averaged 2PCF. In addition, the best-fit value of ^m is 
in quite good agreement with an independent analysis of 
the LRGs by Percival ct al. (2007b) using a power spec- 
trum that yields fim = 0.22 ± 0.04 with a dynamical 
range of 0.01 < k < 0.06/iMpc~^. Our result also agrees 
with tho se from the CMB angular power spectra from 
WMAP CSoer gel et aLllMil . 

We note that the best-fit value of fim is smaller (but 
not significantly so) than the result of E05, fim — 
0.273 ± 0.025. This is not surprising because the ana- 
lytical methods differ from each other. The analysis in 
E05 is based on the spherically averaged 2PCF and uses 
information from smaller scales (s < 40 h~^ Mpc). They 
fix the baryon density parameter as n^h^ — 0.024. We 
obtain a constraint of flmh^ = 0.123lg Q32 from our data 
when the parameter U^h^ — 0.024 is fixed, as listed in 
Table [1] (see also Fig. H]). This result is consistent with 
the E05 result, n,„h'^ = 0.130 ± 0.011. 

All the constraints obtained above are more conserva- 
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Fig. 7. — Likelihood contours for {Qm, f^b, h,ns), assuming a flat ACDM universe. The diagonal panels represent the likelihood functions 
for the four individual parameters with all the other parameters being marginalized over. The other panels show two-parameter constraints 
with the other parameters being marginalized, and each ellipse represents the constraint on the parameter space with 68%, 95%, and 
99% from inward. As in Fig. \E\ the solid (red) and dashed (black) contours are for 40 < s < 200 h~^ Mpc and 60 < s < 160 h~^ Mpc, 
respectively. The best-fit parameters for 40 < s < 200 h~'^ Mpc are Sl„i = 0.218, H;, = 0.0473, h = 0.702, and Us = 1.122 and the minimum 
value of X is Xmin = ^21. 5 with 568 dof. For 60 < s < 160 h'^ Mpc, the best-fit parameters are C™ = 0.208, iit = 0.0462, h = 0.656, 
and ria = 1.030 and x^in = 216.6 with 324 dof. 



TABLE 1 
Summary of constraints on cosmological parameters 





LRG only 
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Note. — The comoving distance s is in units of h ^ Mpc. We assume Q^^h^ = 0.024, w = — 1, 
and rts = 0.98 when they are fixed. 
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Fig. 8. — Likelihood function for a parameter Qrnh^, fitted for 
the scale range of 40 < s < 200 h~^ Mpc. The parameters h, erg, 
and b are marginalized over, and di,h^ = 0.024 and n^ = 0.98 is 



fixed. For the best-fit parameter, x^ 



: 422.2 with 570 dof. 



tive than those in other work using the baryon acoustic 
oscillations because we neglect the small-scale data and 
all CMB information. For example, we obtain worse con- 
straints on h than those on flm- A possible reason for 
this is that discarding the small-scale data makes the 
constraints on flmh^ degenerate along the direction of 
constant acoustic scales (see Fig. 8 of EOS) and this di- 
rection is fairly parallel to the lines of constant O^ but 
intersects with the lines of constant h (Fig. 10 of E05). 

5.4. Dark energy constraint 

Constraining the dark energy is one of the most inter- 
esting applications of the anisotropic 2PCF. As described 
in § |4l we use the information from not only the overall 
shape of the 2PCF but also the geometrical distortion 
with bar yon ridges to constrain the dark energy compo- 
nent. As lMatsuljara fc Szalavl (|2002D indicated, the LRG 
sample is one of the best samples for probing the feature 
of dark energy among current redshift surveys. However, 
it is still difficult for the relatively low-z survey to con- 
strain not only the value of w but also its evolution. In 
this subsection, we therefore assume w to be a constant 
not necessarily equal to —1 and fix flbh^ — 0.024 because 
the baryon density is highly constrained by the analysis 
of the WMAP data fiSpergel et al 2007 ) and big bang 
nucleosynthesis (jBurles et al.ll20QJl ). In Figure [5] we plot 
the joint likelihood functions of J^de and w and obtain 
Ode = 0.770to:o« and w = -O.gstags for the fit to 
40 < s < 200 h~^ Mpc. AU the other parameters but 
flbh^ are marginalized over. Our constraints on dark en- 
ergy parameters are listed in Table [TJ 

We also plot the likelihood function in the bottom left 
panel of Figure B] wher e the no-wiggle power spectrum 
(jEisenstein fc Hij|19 98') is used to calculate the analyti- 
cal formula for the anisotropic 2PCF. Because the dark 
energy parameter is constrained only from the overall 
shape of the 2PCF in this way and degenerates with the 
other parameters without information from the acoustic 
scale, we obtain a poorer fit to the data. Therefore, the 
overall shape is not a dominant effect in constraining the 
dark energy parameter, and the baryon ridges contribute 
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Fig. 9. — Dark energy constraints on Cde and w under the 
assumption of a flat universe and Qi,h^ = 0.024. As in Fig. [7] the 
bottom left panel represents a two-parameter constraint, Q^^-w, 
with h, ris, erg and b marginalized over, and each ellipse represents 
the constraints on the parameter space with 68%, 95% and 99% 
from the inside. The solid (red) and dashed (black) contours are for 
40 < s < 200 and 60 < s < 160 h'^ Mpc, respectively. The two 
diagonal panels represent likelihood functions with all the other 
four parameters marginalized over. The best-fit parameters are 
floE = 0.770, w = -0.93, where x^in = ^21. 2 with 568 dof for 
40 < s < 200 h-^ Mpc, while Cde = 0.786, w = -1.07 and 
X^jj^ = 216.4 with 324 dof. for 60 < s < 160 h~^ Mpc. The 
dotted (magenta) contour in the bottom left panel shows the 68% 
CL when we smooth the oscillatory part of the transfer function 
and use it for the fit of 40 < s < 200 h'^ Mpc. 

as well. 

5.5. Combining with the WMAP results 

So far we have focused on parameter constraints using 
the LRG data only, which are very useful to check its 
result independently, while the obtained constraints are 
inevitably weaker than those we would obtain when other 
data sets are combined. In this subsection we consider 
the additional constr aints using the CMB prior from the 
3 year WMAP data (|Spergel et al.ll2007D . 

We consider two Markov chain Monte Carlo results 
of the WMAP data, w = —1 and constant w cosmolo- 



01 tn e WMAr aata, w = —l and constant w cos 
gies (jTegmarket al.ll2006D . We find 0„ = 0.240^ 



4-1-0.0023 



Ob = 0.04141J; 

^ _ n 7qc+0.050 



h0.023 



0024' " ^ 0-'l°-0.020' 



= 0.947 



0.025' 

-1-0.016 

-0.015' 



and b — 2.15^036 from the former 



0.772to;o33 and w 



0.9715]:^ from 



chain, while Ode 

the latter. These constraints are also summarized in Ta- 
ble [1] and are in very good agreement with the previous 
studies for the joint constraints of the WMAP observa- 
tion with the large- scale structure (e.g., iTegmark et al.l 
120061: iSpergel et aLli2007.) . Although a pure LRG anal- 
ysis cannot tightly constrain O^/i^ or h because of the 
limited range of separations, they are significantly im- 
proved by the prior on the CMB acoustic scale. 

6. CONCLUSIONS 

We have presented the 2PCF in redshift-space for 
the SDSS Luminous Red Galaxy sample considering the 
anisotropy in 2D redshift space. In particular, we have 
focused on the distorted features of the 2PCF in rcdshift- 
spacc from both peculiar velocities of galaxies and geo- 
metrical effect. The distorted features of the Kaiser and 
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finger-of-God effects were clearly detected. The baryon 
ridges, which are the baryonic acoustic features in the 
anisotropic 2PCF, are a nearly spherical object in co- 
moving space. We found indications of baryon ridges 
in the measured 2PCF. Beyond qualitative comparison 
between data and theory, evaluation of the covariance 
matrix is needed for cosmological parameter estimation. 
We constructed the matrix by generating mock samples 
using the second-order Lagrangian perturbation theory 
with an artificial biasing scheme. We have constrained 
the cosmological parameters by comparing the observed 
2PCF with linear theory. 

We have obtained constraints on fundamental cos- 
mological parameters, fi,,,, = 0.218l;nn^y, ^h 
h 



0473+°°^" 



"■'"^-0.117 



and n. 



1 199+0.152 

when we have used the data of 40 < s < 200 h~^ Mpc. 
The constraint on Qm was better mainly because of the 
clear detection of the Kaiser effect, which directly de- 
pends on rim through /3. We have also obtained the 



constraints on the dark energy as floE — 0.770 



+0.051 
-0.040 



and 



w = -O.gataas when we fix f^^/i^ = 0.024 and the other 
parameters, h, Ug, erg, and b are marginalized over. These 
constraints are mainly due to the overall shape of the 
anisotropic 2PCF and the information from geometrical 
distortions including the scale of the baryon ridge. We 
have demonstrated that a pure LRG analysis can con- 
strain w by considering the anisotropy of the structure 
accurately. As for the parameters related to the clus- 



tering amplitude, we have obtained as = 0.66^ 



1.55 



+ 1.42 
-0.75- 



-0.289 
-0.216 



and 



While these two parameters are strongly 



coupled, the degeneracy was alleviated from anisotropy 
of the 2PCF through the redshift distortion factor, /3. 
In addition, stronger constraints on the cosmological pa- 
rameters above were obtained by the CMB prior from the 
3 year WMAP results. All the constraints summarized 
above agree with the previous studies in literature. 

The current analysis can be improved by considering 
two issues below. The first issue is theoretical improve- 
ment accounting for the nonlinearity of the gravitational 
evolution, the redshift distortions and the galaxy bias- 
ing. Although the baryonic signature emerges on very 
large scales, the width of the baryon peak in the 2PCF 
is an order of 10 Mpc. The nonlinearities nontrivially af- 
fect such a feature. In fact, such effects have already 
been investigated using TV-body simulations, higher- 
order perturbation theorie s, and ren ormalizatio n per- 
turba ti on theory (Meiksin et aLN lQQQ: Se o fc EisensteinI 



2005t ISorineel et all l 2005t Ueong & Komats 

Crocce fc Scoccimarrol 120071 : iNishimichi et al.l 

Matsubarall2007f ). The degradation of the acoustic sig- 



nature was well modeled; it was shown that the acoustic 
peak of the linear density fie ld in the 2PCF can be re- 
constructed (jEisenstein et al.ll20 07a,bV In addition, the 
overall shape of the re dshift-space 2PCF is also affected 
by nonlinear dynamics (IScoccimarroll2004 ) . According to 
his result, the redshift-space 2PCF for pairs parallel to 
the line of sight in a random Gaussian field deviates from 
the prediction of standard linear theory even on fairly 
large scales. In this work we do not use the data along 
the line of sight; however, we shall include these issues in 
the future analysis. The theoretical and numerical stud- 
ies also suggest that the biasing i s potentially scale de- 
pendent even on large scales fe.g.. lSchulz &: Whitdl200d : 



ISmith et al.]|2007HColes fc Erdogdull2007f) . which poses a 
serious problem for e stimating cosmological parameters 
from galax v survevs ([Blanton et al.ll2006t iPercival et al.l 
HoOTb; San chez fc Coldl2008D . 

The second issue is the calculation of the covariance 
matrix for the measured 2PCF. As described in § 13.21 we 
have constructed a covariance matrix by the second-order 
perturbation theory. Using N-body simulations which 
fully include nonlinearity provides better estimation of 
the covariance matrix. However, this approach is too 
computationally expensive to produce a large number 
of independent realizations. The approximation by the 
second-order perturbation theory is valid on scales which 
we consider in this work. However, in order to utilize 
the information at smaller scales for more accurate cos- 
mological parameter estimation, we must estimate the 
nonlinearities more accurately using iV-body simulations 
or more sophisticated methods such as a halo occupa- 
tion model from the second-order Lagrangia n perturba- 
tion theory fe.g. JScoccimarro fc Shethll2002D . 

The most important point of our analysis is that we 
directly include anisotropics of the structure. The bary- 
onic features enable to divide the effect of the redshift 
distortions into dynamical and geometrical components. 
The anisotropy due to the geometric distortion, in par- 
ticular, contributes to better estimation of the equation- 
of-state parameter for the dark energy. Various methods 
using the scale of the oscillations as a standard ruler have 
been c onsidered for both the power spectrum and the 
2PCF fEiscnstcin fc Hu'iggS; 'Blake fc Glazcbrook"2003j: 
Hu fc Haiman 2003; Seo fc Eisc iistcin 2003; Matsubar^ 
2004t iGlazebrook fc Blakdl2005f ). This work is the first 



application of the anisotropy in the 2PCF with baryon 
acoustic features to o bserv ational data, which was pro- 
posed bv ' Matsubaral (|2004f) . Direct measurement of the 
growth function from the Kaiser's effect and the two- 
dimensional acoustic scales which depend on Da{z) and 
H{z) is also an attractive challenge of the analysis using 
the anisotropic 2PCF, so these topics will be definitely 
pursued in future work with an improved LRG sample. 
The baryonic signature from the redshift range of SDSS 
LRGs is not strong, because their number density is rel- 
atively small and nonlinear effects weaken the baryonic 
feature. There are many plans for constraining the dark 
energy by future wide- field, deep galaxy survey s: the 
Fiber Multiobject Spectrograph (FMOS; Kimura et al.l 



2003h. Wide-Field Multiob.ject Spectrogr aph (WFMOS; 
Glazebrook fc Blakd l2005t Bassett et al.l[2005.). B arvon 



Oscillation Probe (BOP; IGlazebrook et al.l 120051 ) . and 
the Hobby-Ebe rly Dark Energy Experiment (HETDEX; 
iHill et aLll2004f ). and so on. When the baryonic signa- 
ture is detected with high accuracy from future redshift 
surveys, the analysis of the anisotropic 2PCF as in this 
work will be an important ingredient for stringently con- 
straining properties of the dark energy. 
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